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Abstract. We introduce and study the toric fiber product of two ideals 
in polynomial rings that are homogeneous with respect to the same 
multigrading. Under the assumption that the set of degrees of the vari- 
ables form a linearly independent set, we can explicitly describe gener- 
ating sets and Grobner bases for these ideals. This allows us to unify 
and generalize some results in algebraic statistics. 

1. Introduction 

A common problem in algebraic statistics is to convert the parametric 
representation of a statistical model into the implicit representation in terms 
of finding the defining prime ideal of the model. This is a special case of 
the implicitization problem that arises frequently in computational algebraic 
geometry. In algebraic statistics, we are usually presented with a family of 
statistical models and we would like to find a theorem which gives a complete 
description of all the ideals for all the statistical models in this family. A 
useful approach has been to try to find decomposition rules for the models 
and the resulting ideals, and subsequently reduce the problem to finding 
the defining prime ideals in a few special cases which can then be handled 
theoretically or using a computer algebra system. This approach has played 
a role in attacking the problem of determining phylogenetic invariants for 
various tree-based models of evolution ^ EI and for studying Markov bases 
of hierarchical models |H1 E| ■ 

In this paper, we introduce the toric fiber product, an operation that 
takes two homogeneous ideals with compatible multigradings and produces 
a new homogeneous ideal. This operation generalizes the Segre product of 
two schemes as well as the gluing operations for toric ideals that appear in 
[HI El E|' When the underlying grading group has special structure, we 
are able to explicitly compute generating sets and special Grobner bases 
for the toric fiber product from generating sets and Grobner bases of the 
component ideals. 

Given a positive integer n, let [n] = {1,2, ... ,n} denote the set of the 
first n positive integers. Let r > be a positive integer and s,t G Z'^q be 
two vectors of positive integers. Let 

K[x]=K[xi \iG[r],j e[s^]] 

and 

K[y]=K[yi \ze[r],keM 
1 
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be multigraded polynomial rings subject to the multigrading 

deg(4) = deg(y|) = G Z'^. 

We assume throughout that there exists a vector lo & Q'^ such that w-^a* = 1 
for all i. This implies that ideals in K[x] or K[y] that are homogeneous with 
respect to the multigrading are homogeneous in the usual sense. Denote by 
A = {a\ . . . , a'"} and let NA be the affine semigroup generated by A. 

If / and J are homogeneous ideals (with respect to the multigrading) in 
K[x] and respectively, the quotient rings R = K[x]/I and S = K[y]/J 
are also multigraded rings. Let 

K[z]=K[zi,, \i€[r],j e[si],keM 

and let 4>i^j : K[z] R (8)k S be the ring homomorphism such that zjj. 

x'j yl- 

Definition 1.1. The toric fiber product of I and J, denoted / J is the 
kernel of (pi^j: 

Ixj^J = kBY{(t):^j). 

Two fundamental examples, illustrating the coarsest and finest possible 
multigradings, are Segre products and sums of monomial ideals. 

Example 1.2. Suppose that r = 1 and / = J = 0. Then (f)j^j is the ring 
homomorphism 

(t>i,j : K[zjk \je[s],ke [t]] ^ K[xj,yk \je[s],ke [t]] 

Zjk ^ Xjyk. 

The toric fiber product is 

I x_aJ = ker(0/,j) = {zj^k-^Zj^ki - Zj^ki^jiki I jij2 G [s], ki, k2 G [t]); 

that is, / x_4 J is the ideal of 2 x 2 minors of a generic matrix. 

Example 1.3. Suppose that r > and s,t = are the all ones 

vector. Suppose that deg(x*) = deg(y*) = deg(^*) = Cj, the ith standard 
unit vector. If / G and J € are homogeneous with respect to 

this multigrading, they must both be monomial ideals. Then the toric fiber 
product is simply 

IxaJ = I{z) + J{z) 
where I{z) denotes the ideal / with z variables substituted for the x variables 
(and similarly for J{z)). 

Our main interest in toric fiber products is when I and J are the prime 
ideals of unirational varieties. Then the ideal / x_4 J is also prime and 
defines a unirational variety. In practice, we are often presented with the 
parametrization of a unirational variety and we are interested in finding its 
defining ideal. One useful tool is to find a nice grading such that the ideal 
is, in fact, a toric fiber product. This grading is usually considerably coarser 
than the finest grading associated to the ideal. If .4 is a linearly independent 
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set, we can determine generators and Grobner bases of the toric fiber product 
/ x_4 J explicitly from generators and Grobner bases of I and J. 

Example 1.4. Let (j) be the ring homomorphism 

: ]K[gij,j3,jg I ii,Z3,i5 G [3]] ¥\ai^i^,hi^i^,Ci^i^,di^i^ \ 11,13,15 G [3],i2,iA ^ 

2 2 

12 = 1 14 = 1 

Each of the polynomials appearing in the parametrization is homogeneous 
with respect to the grading deg(pij,i3,jg) = 6^3, the ia-th standard unit 
vector. The polynomials appearing in the parametrization can be written 
in factored form as 

\i2 = l / \j4 = l / 

The ideal K = ker((/)) is a toric fiber product K = I I ■ The underlying 
ring R = K[qi-^^,i^ \ 11,13 G [3]] has the grading deg(gj^j3) = 6^3, the kth 
standard unit vector. The ideal / is the kernel of the ring homomorphism 

K[ai^i2,bi2i^ \ 11,13 & [3],i2 G [2]] 
2 

j2=l 

Thus, I is the principal ideal generated by the determinant of the matrix: 

Qui QU2 QU3 

'?2»l 92.2 Q2»3 
K^Q3»1 Q3»2 g3»3y 

Using the machinery in Section [3 one can show that the ideal K = I Xj^I 
is generated by determinants of flattenings and slices of the 3-dimensional 
tensor (qiimismis)- In particular, K is generated by the 3x3 minors of the 
matrices 

QUUl QUU2 QUU3 QU2»1 QU2»2 QU2»3 9l»3»l QU3»2 9l«3»3 ' 
Q2»U1 Q2»U2 Q2»U3 Q2»2»l Q2»2»2 Q2»2»3 Q2»3»l Q2»3»2 Q2»3»3 
yQStUl Q3tU2 Q3tU3 Q3t2»l Q3t2»2 93»2»3 93»3»1 93»3»2 93»3»3/ 

ft.l.l '?1.2»1 92.1.1 '?2.2»1 92»3»1 93»1»1 Q3»2»l Q3»3»l 

9l.l»2 '?1.2»2 '?1.3»2 Q'2.1»2 '?2.2»2 92»3»2 Q3»U2 93»2»2 93»3»2 
^Ql^l^S 91.2.3 QU3»3 Q2»U3 Q2»2»3 Q2»3»3 Q3»U3 Q3»2»3 93»3»3y 

together with the 2x2 minors of the matrices 

QUUl QUU2 QUU3\ I QU2»1 qU2»2 qU2»3 ' 
92»1»1 92»1»2 92»1»3 j I Q2»2»l Q2»2»2 Q2»2»3 
i93»l»l Q3»U2 93»1»3/ \'/3»2»l 'Z3»2»2 '?3»2»3/ 
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9l»3»l 9l»3»2 9l.3»3 \ 
92«3»1 '?2»3»2 '?2.3»3 j • 
'?3«3»1 '?3.3»2 '?3.3»3 / 

Furthermore, this collection of 2 x 2 and 3x3 minors form a Grobner basis 
for K. This example is a special case of Corollary 13.81 □ 

The main focus of this paper is on the special case of toric fiber products 
where ^ is a linearly independent set. As we will see, this played a significant 
role in Example 11.41 It is probably impossible to recover explicitly the 
generating set of -fC = / x_4 J from the ideals I and J, if A is not linearly 
independent. Indeed, in the next section, we will see an example where I is 
generated by quadrics and J = (0) but I J requires minimal generators 
of arbitrarily large degree. 

The outline for this paper is as follows. In the next section, we use 
a contraction of an ideal under a monomial homomorphism to determine 
generating sets and Grobner bases for the toric fiber products I x^J, when 
A is linearly independent. In Section |31 we consider applications of the 
main result. This includes proofs of some known results about the defining 
ideals of products of projective schemes, and their Grobner bases. We also 
illustrate how the toric fiber product arises in algebraic statistics. This 
allows us to unify results about reducible hierarchical models and group- 
based models on phylogenetic trees. 

2. Contractions Under Monomial Homomorphisms 
Let B € Z"^^" be a dxn integral matrix. Consider the ring homomorphism 

CpB ■■ ^[Zl, ...,Zn]—^ K[ti, ...,td\ 
i=l 

We call (j)B a monomial homomorphism. The ideal Ib = ker((/)B) is called a 
toric ideal. Sturmfels' book jl6j is a standard reference for background on 
toric ideals. In this section, we consider the contractions of arbitrary 

ideals I in and apply these results to toric fiber products. The main 
idea here is to compare the initial ideals of I to the initial ideals of cj)~j^{I). 

Let uj E Z>g be a vector of weights. The vector uj induces a partial order 
on the set of monomials in by declaring that -<^^ if uj'^a < uj'^h. 
The partial order ~<i^ is called a weight order on M.[t]. Note that it need not 
be a term order on 

Given a polynomial / G the initial form in^(/) is the sum of all 

terms of / that have the highest weight with respect to the partial order 
-<(^. If I is an ideal of the initial ideal in^(/) is the ideal 

in^(/) = (in^(/) I / G /) . 

Our main use for for weight orders comes from the following useful fact. 
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Proposition 2.1. jl61 Prop. 1.11] For any term order -< and any ideal 
I C ]K[t] there exists a vector uj € Z>q such that m.i^{I) = in^(/). 

We say that a finite collection of polynomials G C / is a Grobner basis 
of I with respect to the weight order oj if (in^ {g)\g ^ G) = in^; (/) and this 
ideal is a monomial ideal. A collection of polynomials such that (in^((7) | g G 
G) = in^(I) is called a pseudo-Grobner basis. Any Grobner basis or pseudo- 
Grobner basis of / generates /. 

Every weight order on K[t] determines a weight order ^(f)*^uj on K[z] 
via the pullback of lu through (p. That is, (j)*^io = co'^B. Note that this 
construction has two important properties. First, if is a monomial then 
its weight with respect to ^^u; equals the weight of 4>b{z^) with respect to 
UJ. Second, if / G ker(0B) is a binomial then in^|^(^(/) = /, which implies 
that va.fj,*^^{lB) = Ib for all u; G R'^. 

Lemma 2.2. Let I be an ideal in ]K.[t]. Then 

in^^^(</.^i(/))C</,^i(in^(/)). 

Proof. Let / G (/)^^(/). We must show that (t>B{^^<i)*gUj{f)) G Without 
loss of generality, we may assume that / is reduced with respect to any 
Grobner basis of the toric ideal Ib- This is because Ib C in^* for 
any / and ^b{Ib) = 0. Since / is reduced with respect to Ib, this means 
there is at most one term of / in each B graded degree. In particular, 
there can be no cancellation amongst the terms of (f)B{f). Since each pair 
of monomial of / and image monomial of (pBif) have the same weight with 
respect to (p*^uj and m respectively, we deduce that (pBi^^rjy'gOjif)) ^ ina;(/) 
which completes the proof. □ 

Lemma 2.3. Let M = {mi, . . . ,mr) C be a monomial ideal. Then 

cP^HM) = ^b'H^i)) + ■■■ + ^pB^i'^r)). 

Furthermore (f)^^{M) = M' + Ib where M' is a monomial ideal. 

Proof. We use the same argument as the proof of Lemma 12.21 In particular, 
suppose / G (f)~j^{M). We may suppose that / is reduced with respect to 
any Grobner basis of Ib- Then each monomial in / maps to a monomial 
in IfC[t], and / G (f)~^{M) if and only if each monomial of (psif) belongs 
to M. This means that if n is a monomial of /, (j)B{'n^) G ("^j) for some i. 
Furthermore, this shows that each monomial of / belongs to (t)~j^{M). □ 

Lemmas 12 . 21 and l2 . 31 suggest a strategy for determining the ideals (j)'^{I). 
First, we compute an initial ideal of I. Then we determine cf)'^ {in^^{I)) 
using combinatorial arguments. Then, if we are lucky, we find a collection 

of polynomials G C <j)~^^{I) such that ^in0^^(G)^ = m^*^i^{(j)']^^{I)). Then 

we can conclude that G is a pseudo-Grobner basis for cl)^{I) with respect 
to the weight order -<4,*gU)- 
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In general, this strategy is not possible to implement, because either 
m^*^^{4>~^^ {!)) 7^ ((>]^^{m^{I)) or there is no combinatorial description of 
4>^{M) where M is a monomial ideal, or both. However, in the special 
case that arises when taking a toric fiber product, we will see that there is 
a simple answer to both problems. 

To give our main algebraic result concerning the Grobner bases and gen- 
erating sets of the toric fiber products / x_4 J, we first need to show that 
the toric fiber product is the contraction of a monomial homomorphism 
<i>~Q'{I + J), for suitable B. To this end, we will derive an alternate descrip- 
tion of the toric fiber product which fits into this framework. 

With the setup from Section^ consider the monomial homomorphism 



(1) 4fc ^ ^^fc- 

For the remainder of this section B denotes the matrix arising from a 
toric fiber product according to Equation ^ This matrix only depends on 
r, s, and t. 

Proposition 2.4. 

/X^ J = </.-!(/ + J). 

Note that the ideal I + J is considered as an ideal in y] after taking 
the extensions of / and J. 

Proof. Note that IfC[x, y] ^ K[x] (g)KlfC[y]. Hence + J) ^K[x]//®k 

]K[y]/J. Given any ring homomorphism (j) : R ^ S and induced ring homo- 
morphism (p : R ^ S/I, we have ker((/)) = (p~^{I), and this completes the 
proof. □ 

Lemma 2.5. Let m = x^}x^? ■ ■ ■ be a monomial in y]. Then 
rB\{rn)) = i^,^^^^ ■ ■ ■ zl,^ I G [t,J, . ..k,e + Ib. 

Similarly, if n = y^^^y^^ ' ' ' ^fc^ ^ monomial in K[x, y], then 

^sHi^)) = {^k,^k, ■ ■ ■ ^uk, I h e ...jd^ K]) + Ib. 

Proof. It suffices to handle the first case. By the second part of Lemma 
12.31 it suffices to determine the monomials that belong to 0^^((m)). Denote 
by M' the monomial ideal on the right hand side of the equation. Given 
a monomial m' = 11-^]; fc; ^ "^[A^ image is the monomial H^j'yfcj which 
belongs to (m) if and only if Y\x^j^ belongs to (m). But this implies that 
these exists a monomial in M' dividing m' . □ 



TORIC FIBER PRODUCTS 



7 



Proposition 2.6. Let B he the matrix representing the monomial homo- 
morphism arising from the toric fiber product. Then 



Ib 



7 7 
il^l ^2^2 



1 < z < r, 1 < ji < j2 < s j , 1 < A;i < ^2 ^ 



and these quadrics are a Grobner basis for Ib with respect to any term order 
that selects the underlined terms as leading terms. 



Denote by Quad^ the set of quadrics described in Proposition 127 

Proof. First of all, there exist term orders which select the underlined terms 
as leading terms. Indeed, let ~< be the lexicographic term order such that 
z]\ki -< if h < h or ii = 12 and ji < j2 or ii = 12 and ji = j2 and 
ki > k2. In particular, -< selects the underlined terms of the quadrics in 
Quad^ as leading terms. 

Since Ib is a toric ideal, it suffices to show that if / is any binomial in Ib, 
there exists a quadric g G Quad^ such that \n^{g)\\n^[f). To each binomial 



/ 



^d 

jdkd 



we associate the tableaux of indices 

ii ii kx 

f 



id 3d 



kri 



fd 



k[ 



kL 



Note that two individual tableau represent the same monomial if and only if 
one can be obtained from the other by swapping rows. A binomial / belongs 
to Ib if and only if the image of each monomial under (/>_b is the same. In 
tableau notation, this can be expressed as 



h 


ji 




k 






'A 3[' 




i[ k[ 


_ id 


3d _ 






kd _ 




Jd j'd. 




7' k' 

.d '^d _ 



In each expression, the first tableau denotes the indices of the x variables 
and the second tableau denotes the indices of the y variables. Thus, after 
rearranging the rows of the second tableau of / we may write 



/ 



k ji h 



I'd 3d 



kri 



k 31 k[ 



I'd 3d 



k' 
'^d 



Furthermore, we will assume that the rows of each tableau are ordered lexi- 
cographically with 111 ~< 112 -<■■■. A monomial is reduced with respect to 
the set Quad^ if and only if there is no sub-tableau of the form 

i ji k2 
i j2 ki 

where j\ < j2 and ki < ^2- This implies that in a monomial that has been 
reduced with respect to Quad^, for each fixed value of i, both the j's and 
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/c's are nondecreasing down their column in the subtableau containing ah 
the rows with it = i. However, among all the tableaux that have the same 
image under (ps, there is only one that has this property and it is minimal 
with respect to Thus, if / G Is its leading term must be divisible by 
some leading term in Quad^. □ 

For the rest of this section, we will assume that v4 is a linearly indepen- 
dent set. Let / G / be a homogeneous polynomial (with respect to the 
multigrading by N.4). If / has degree d and v terms we can write 

V 

J, \ ^ if q 
u=l 

where each G K is a coefficient. However, the fact that / is homogeneous 
and A is linearly independent guarantees that each multiset of upper indices 
Mu = {ii, ■ ■ ■ , i^} is independent of u. That is, for all u, u', M„ = M„/. Thus 
after possibly rearranging the indeterminates appearing in each monomial 
we can always write 

V 

J" \ ^ ii i2 id 
u=l 

Now let = (fci, . . . , kd) with ki G [U^], k2, G [U^], . . . , A;^ G [Uj and consider 
the polynomial fk G K[z] defined by 

V 

u=l 

Note that since f & I, the new homogeneous polynomial fk & I x a J ^'^^^ 
all k. This follows because 

Mfk) = yi\yg---yZ-fei. 

Definition 2.7. Let A be linearly independent and let F C / be a collection 
of homogeneous polynomials. To each / G -F we associate the set Tf = 

nf=i[^ij of indices. Denote by 

Lift(F) = {/fe \feF,keTf} 

which we call the lifting of F to / J. If G C J is a collection of homo- 
geneous polynomials we define Lift(G) in the analogous way. 

Theorem 2.8. Suppose that A is linearly independent. Let F <Z I he a 
homogeneous Grobner basis for I with respect to the weight vector loi and 
let G G J be a homogeneous Grobner basis for J with respect to the weight 
vector L02 ■ Then 

Lift(F) U Lift(G) U Quads 
is a pseudo- Grobner basis for I x ^ J with respect to the weight order ^^(a;i,a;2). 
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Proof. First of all, note that since the generators of / and J are in disjoint 
sets of variables, FuG is a Grobner basis for I+J with respect to the weight 
order {loi,uj2)- Let M = in(^^ ,^2)(-^ + J) be the initial ideal of I + J with 
respect to (cji, 0^2)- Since M is a monomial ideal and each minimal generator 
of 1 + J belongs to ]K[a;] or we can use Lemmas 12 . 31 and 12 . 51 to compute 
0^ (M). However, each of the monomials appearing in 0^ (M) appears 
as the (p'^(uJi,uJ2) leading term of some polynomial in Lift(F) or Lift(G). 
Additionally, Quad^ generates Ib- So we deduce that the initial forms in 
Lift(F) U Lift(G) U Quads C I x^J generate 0b^(M). Thus, by Lemma 
12.21 we deduce that (l)]^^{M) = in^*^(^^ '^A J) and this completes the 
proof. □ 

Theorem 2.9. With the same assumptions as Theorem \2.S\. let uj be a 

weight vector such that Quad^ is a Grobner basis for Ib ■ Then 

Lift(F) U Lift(G) U Quads 

is a Grobner basis for I Xj^^ J with respect to the weight order 0s('^i)'-^2) + e'^ 
for sufficiently small e > 0. 

Proof. Ifwe choose e very smah, we will have in0|^(^j^^2)(/fc) = m.4>*^{uji,uj2)+eMk) 
for all fk G Lift(F) U Lift(G). This implies that 'va.^*g(^^^,ui2)+eui{I J) = 
in^(in^*^(^^ J)) since we only need to determine the initial terms 

of (Quads). But because oj induces a term order that make Quad^ into a 
Grobner basis for Ib we are done. □ 

Corollary 2.10. Let A be a linearly independent set. Let F be a homo- 
geneous generating set for I and G be a homogeneous generating set for J. 
Then 

Lift(F) U Lift(G) U Quads 
is a generating set for I Xj^^J. 

Proof. If F and F' generate /, Lift(F) and Lift(F') generate the same ideal. 
This holds, in particular, if F' is a Grobner basis for /. A similar statement 
holds of G as well. Thus 

(2) (Lift(F) U Lift(G) U Quads) = (Lift(F') U Lift(G') U Quads) 

where F' and G' are Grobner bases for / and J respectively. But the ideal 
on the right hand side of Equation |21 is I Xjs,J , since a Grobner basis of an 
ideal generates that ideal. □ 

Corollary 2.11. Suppose that A is linearly independent, and that m.i^^{I) 
and m.i^^{J) are squarefree monomial ideals. Then in(^|j(aji,a;2)+<:'^(-^ ^A J) *s 
a squarefree monomial ideal. 

Proof. Since 'm^-i^(I) and in^jC"^) squarefree monomial ideals, then the 
(j)*^{uj 1,102) + euj initial term of every polynomial in Lift(F) and Lift(G) is 
squarefree. Since the w leading terms of Quads a-re also squarefree, and the 
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union of Lift(F), Lift(G), and Quad^ form a Grobner basis for / x_4 J we 
deduce that in^*^(^j ,^2)+e<^(-^ J) is a squarefree monomial ideal. □ 

It seems natural to ask whether any more refined information about the 
toric fiber product / x_4 J can be computed from / and J (for instance, 
homological properties, Betti numbers, Hilbert series). In view of Example 
11.31 it seems that there is no hope of explicitly determining any of these 
properties from / and J, since the sums of even simple monomial ideals can 
display complicated behavior. There is, however, a compact description of 
the multigraded Hilbert function and series of I x_4 J. We refer the reader 
to for an introduction to multigraded Hilbert functions. Note that since 
we have assumed that the underlying grading is positive, dimK((IK[x]//)u) 
and dimK((IK[y]/J)„) are finite for all u € NA. Thus the Hilbert function 
and series exist. 

Corollary 2.12. Suppose that A is linearly independent. Then the Hilbert 
functions ofM.[x]/I, K[y]/J and x_4 J satisfy 

h{K[z]/I x_aJ-u) = h{K[x]/r,u)h{K[y]/J;u). 

Thus the Hilbert series ofK.[z]/I x_4 J is the Hadamard product 

H{K[z]/I x^ J;t) = i7(lC[x]//;t)(lC[y]/J;t). 

Proof. It suffices to produce the first equation, since this is the definition of 
the Hadamard product of two series. Let G be a Grobner basis for / x_4 J 
constructed according to Theorem [221 Let M = in^^^^^^ (^2)+eaj (I J) be 
the corresponding initial ideal. We need to count the number of monomials 
in the u graded pieces of \ Af . A monomial m = belongs to ]K[2;] \ M 
if and only if m is not divisible by any of the initial terms of polynomials 
of Lift(-F), Lift(G) or Quad^. This implies that (j)B{m) is not divisible by 
the leading terms of polynomials of F or G. The monomial (j)B{m) is a 
product of a degree u monomial in ]K[x] \ ini^^(/) and a degree u monomial 
in K[y] \ id^^^J). Any such product has a unique preimage that is not 
divisible by any leading term in Quad^. Thus we have shown that the 
standard monomials of / x_4 J of degree u are in bijection with a product 
of a standard monomial of / of degree u and a standard monomial of J of 
multidegree u, which yields the desired equality of Hilbert functions. □ 

To close this section, we provide an example which shows the necessity 
of the condition that ^ be a linearly independent in all of the preceding 
theorems. 

Example 2.13. For a fixed vector of positive integers d = (^1,^2,(^3) with 
di > 1 for all i, consider the ring homomorphism 

</>d : Klpi^i^i^ I ij G [dj]] ]K[ajii2, ftj^ig, ^343 | ij G [dj]] 

Piii2i:i ' ^ '^4142 ^ii 43 '^j2*3 ■ 

Computing the minimal generators of the ideal = ker(0rf) for various 
values of d is a benchmark problem for algorithms for computing Grobner 
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bases of toric ideals ^1 . These ideals are extremely complicated and 
while there are explicit generating sets known for some special values of d, 
there is, at present, no uniform description of the generating sets of these 
ideals. Work of De Loera and Onn suggests that it is impossible for 
any simple description of a generating set to exist for all d. In particular, 
it is known that Kd requires minimal generators of degree at least than 
2min((ii,d2,(i3) 

The ideals are examples of toric fiber products that do not have lin- 
early independent A. In particular, consider the grading deg(pj^i2j3) — 
Gil ® where denotes the k-th standard unit vector. Consider the 
ring homomorphism 

jpd ■■ IKfeiiais h ^ [dj]] IK[6iij3,Q2i3 I ij G [dj]] 

and let 1^ = keT{ipd)- This ideal is generated by quadrics: it is an ideal 
of the form Ib arising from a toric fiber product with linearly independent 
A as in Proposition 12.61 Grade K[q] by deg(qjjj2«3) = ® ^i2- ^et = 
(0) be the zero ideal in lC[o] with the grading deg(ajjj2) = ® ^12- 
A = {cj^ © ej2 I ij € [dj]}- Then K^. = J^- So although 1^ and 

are generating in degree 2 or less, can require minimal generators of 
arbitrarily large degree. 

3. Applications 

3.1. Segre products. The simplest example of the toric fiber product is 
the usual Segre product of two projective schemes. In this setting I C IfC[x] 
and J C ]K[y] are homogeneous ideals in the usual coarse grading by degree. 
The monomial homomorphism 0b is 

^jk ' ^ -^iUj ■ 

In this case, all variables Xj and yk have degree 1 and since ^ = {1} is a 
linearly independent set we are in a position to apply the results of Section|2l 
Theorem 12 . 81 shows how to produce Grobner bases for / J from Grobner 
bases for / and J. For instance, we deduce the following corollary, which 
appears in [TT] . 

Corollary 3.1. Let d be a vector of nonnegative integers greater than 1. Let 
Kfi he the defining ideal of the product of projective spaces P'^i^^ x • • • x P'^"^! 
under the standard Segre embedding into p'^i'"'^"^!. Then Kd is generated 
by the 2x2 minors of all flattenings of a generic di x • • • x d„ tensor. 
Furthermore, these 2x2 minors form a Grobner basis for and the ring 
W^z\/Kd is Cohen- Macaulay. 
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A flattening is a matrix obtained by partitioning the n indices of a tensor 
into two nonempty sets. For instance, a flattening ofa2x2x2x2 tensor 
with partition {1,2}|{3,4} is 

'^^llll ^1112 ^1121 ^1122\ 

^1211 ^1212 2;i221 ^1222 

^2111 ^2112 2:2121 ^2122 

\^2211 ^2212 2:2221 2:2222/ 

In ^] , Ha refers to the 2x2 minors of flattenings as minors of a box-shaped 
matrix. 

Proof. It is easy to see that any 2x2 minor of a flattening belongs to 
Kii and, furthermore, these are the only degree two binomials in the ideal 
Kd- Note that = K'^ x^ (0) where d! = {di, . . . ,dn-i) and A = {!}. 
By repeatedly applying Theorem 12.81 we can construct generating sets and 
Grobner bases for Kd. The lifting operation preserves degrees and every 
polynomial in Quad^ has degree two, so the resulting Grobner bases have 
degree two. Since the only degree two binomials in are the 2x2 minors 
of flattenings, these must form a Grobner basis. By Corollary 12.111 the 
resulting initial ideal is squarefree. Since Kd is a toric ideal, this implies that 
the K.[z]/m^{Kd) is Cohen-Macaulay (the simplicial complex associated to 
in^(/Crf) is a regular triangulation 16, Theorem 8.3] and, hence, a shellable 
ball) which in turn implies that ]K[2:]/i^(i is Cohen-Macaulay. □ 

3.2. Reducible models. Probably the first instance where Theorem 12.81 
was used in some generality was in the study of the class of reducible hier- 
archical models in |S] and Hierarchical log-linear models are a class of 
statistical models use in the analysis of multivariate discrete data. To each 
such hierarchical model is associated a toric ideal lA,d- The generators of 
the toric ideal /a,^ are useful for performing various statistical tests, as first 
demonstrated in [7j. In this section, we will only describe these models in a 
purely algebraic language and show how results about the Grobner bases of 
reducible models follow from the theory in Section [21 

Let A be a simplicial complex with ground set [n] and let d = (di, . . . , dn) 
be a vector of integers with di > 2 for all i. We suppose that |A| = Upfz^F = 
[n] . For a subset F C [n], we use the notation Dp to denote the set of indices: 

DF=l[[dk]. 

For a given string of indices i G -D[„], If is the subvector ip = (iki, ■ ■ ■ ,iks) 
where F = {ki, . . . , kg}. Denote by 

K\p] =K[pi\ie 

and 

K[a] = K[a[^ \ F =G facet(A),jF G Dp]. 
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Consider the ring homomorphism 

n 

Fefacet(A) 

Definition 3.2. The toric ideal /^.d = ^^^i4'A,d) is the ideal of the hierar- 
chical model defined by A and d. 

Example 3.3. Let A be the simplicial complex with facets facet(A) = 
{{1,2}, {2, 3}, {3, 4}}. Then 

0A,d : Klpi^i^i^i^ : ii G [di]] > ojajf I Jl ^ [dl]] 

{1,2} {2,3} {3,4} 
Kl«2i3*4 ' ^ '^1112 ^12*3 ^^1314 " 

If A is the simplicial complex with facets facet(A) = {{1, 2}, {1, 3}, {2, 3}}, 
then we get the ring homomorphism from Example 12.131 □ 

Definition 3.4. A simplical complex A is called reducible if there are two 
subcomplexes Ai, A2 C A such that Ai U A2 = A and Ai n A2 = 2"^ for 
some S C [n]. The set S is called a separator. 

For instance, the first simplicial complex from Example 13.31 is reducible 
with S = {1} or 5 = {2}, whereas the second simplicial complex from 
Example 13.31 is not reducible. We will show that if A is reducible, the ideal 
lA,d can be written as a toric fiber product. This will allow us to deduce 
Theorem 4.17 from |13j . 

To this end, let d^ and d^ be the induced subvectors of d on the index sets 
|Ai| and IA2I respectively. That is d^ = (i|Ai| and = d\^^\. For / G {1,2} 
let 

K[p\i = IfC[pi|^^| I G Z)|A,|] 
and consider the ring homomorphism 

Fefacet(Ai) 

We denote by /a;,^ the kernel of (pAi^h which is the toric ideal of the 
hierarchical model associated to A;. Since Ai n A2 = 2'^, there exists two 
facets, Fi G Ai and F2 G A2 such that S C Fi and 5 C F2. We introduce 
a grading on K[a] so that 

deefa^ ) - / ''O^^s if ^ ^ {^1, ^2} and 
^«S^"jF^ I otherwise 

The vector Ci^ for G ds is the standard unit vector in Z^" with a 1 in 
the position and a zero elsewhere. This multigrading on K[a] induces 
a multigrading on lC[p], and IK[p]2, where, for instance, we take the 

degree of pi to be deg(pi) = Cj^. Thus, all of lA,d, -^Ai,di and 1/^24'^ 
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homogeneous with respect to this multigrading, because the maps (pA,d, 
0Ai,di s-iid 4>i\2,d? ^-ll preserve the multidegree. 

Theorem 3.5. Let A he reducible with components Ai and A2 and separator 
S. Then 

with A linearly independent. 

Proof. It is clear that A is hnearly independent since it is a coUection of 
disjoint standard unit vectors. Suppose that T is an arbitrary face of A. 
We first note the general fact that if we modify the ring homomorphism (/>A,d 
so that 

Pi^O'iT- n "'ip 

Fgfacot(A) 

this does not change ker((/)A,d), since any variable af^ must appear with 
precisely same multiplicitly as for any facet F with T C F in the image 
of a monomial (j)A,d- Consider the modified parametrization 

Pi -(4)^- n < 

Fefacet(A) 

where S is the separator. We factorize the expression on the right as: 

K)'- n < = n U- n < 

FGfacet(A) ^£{1,2} \ FGfacet(A)nfacet(A() 

Since the expression on the right inside the parentheses involves a product 
of terms for all facets of A^, this product represents the parametrization for 
^Ai,d'- other words, we have shown that the ring homomorphism (j)/\^d 
factors through 

Thus /A,d is a toric fiber product. □ 

3.3. Reducible models with hidden variables. Pushing the idea from 
Section 13.21 one step further, we can also use the machinery to compute 
the ideals of reducible models from submodels when some of the random 
variables are hidden. Parametrically, we have the following setup for the 
algebraic description of a hidden variable models. 

Let A be a simplicial complex on [n], d = {di,...,dn) and vector of 
integers with di > 2 for all i. Let H = {hi, . . . ,ht} C [n] be the collection 
of hidden nodes and O = [n]\H he the collection of observed nodes. Let 

K[g]=lC[gi I ioeDo,ii = »ifleH] 
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and let K\p], (j)A,d, and 1^4 be defined as in Section ESI If io G and 
in £ Dh we use the notation Pi^j^ to denote the indeterminate pi such that 
il = (io)/ if / G O and i; = (j_f/)i if I £ H. Consider the ring homomorphism 

Denote by IhA4 = ker('i/'_H' o which is the ideal of the hidden vari- 

able hierarchical model. This ideal is rarely a toric ideal. Hidden variable 
graphical models have been studied in jH] and ^HI from the perspective of 
computational algebra though we seem to be the first to write down some 
general principles for determining their defining prime ideals. 

Definition 3.6. We call the hidden variable ideal IhaA reducible if A is a 
reducible simplicial complex and H D S = f/>. 

Suppose that lH,A,d is a reducible hidden variable ideal. Let Ai and A2 
be the two component subcomplexes, let Hi = H n |Ai| and H2 = H n IA2I 
and let and be the two induced vectors of indices. Denote by 
and K[q]2 the two rings with variables indexed by the elements of and 

Theorem 3.7. Let lH,A,d be a reducible hidden variable ideal. Let S be the 
separator. In each of the rings and ^[q]2, let the degree of a 

variable be deg(gi) = Cig. Then 

with A linearly independent. 

Proof. The same proof as of Theorem 13.51 applies here. □ 

As an example of an application of Theorem 13.71 we will deduce a Corol- 
lary generalizing Example 11.41 We call this example the partially hidden 
Markov chain. Let A = {{1, 2}, {2, 3}, {3, 4}, . . . , {2n, 2n + 1}} be a chain 
of odd length and suppose that H = {2, 4, 6, ... , 2n} consists of all the even 
numbers. To describe the generators of the ideal Ih,A4 we need two matrix 
constructions. First, each even number 2j G [2n -|- 1] defines a flattening 
of the n + 1 tensor (gj) into a matrix Xj. The rows and column indices of 
the matrix Xj are the elements oi Dp with i^ = {l,3,...,2j — 1} and the 
elements of Dp' with F' = {2j -|- 1, 2j -f- 3, . . . , 2n -H 1}, respectively. Thus 
Xj is a did^ ■ ■ ■ c^2j-i x '^2^+1^2^+3 • • • c^2n+i matrix. The entry in the ip row 
and jpi column is Qipjp, (with appropriate w's added). 

Second, to each odd number 2j + 1 G [2n -|- 1] with < j < n, and 
each i G [ci2j+i] we introduce a matrix Yj i which is a flattening of an n 
dimensional slice of the (n-l-l)-dimensional tensor {q{). The row and column 
indices of the matrix Yj^i are the elements of Dp with F = {1,3,... , 2j — 1} 
and the elements of Dp' with F' = {2j + 3,2j + 5, . . . ,2n + 1}. Thus, Yj,i 
is a did3 ■ ■ ■ d2j-i x ^2^+3^2^+5 • • • '^2n+i matrix. The entry in the ip row 
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and jp/ column is qipijp, (with appropriate bullets added. Examples of Xj 
and Yj^i are illustrated in Example 11.41 (note that the second 3x9 matrix 
in Example 11.41 is the transpose of X^). 

Corollary 3.8. Let A = {{1, 2}, {2, 3}, {3, 4}, . . . , {2n, 2n + 1}} be a chain 
of odd length and suppose that H = {2, 4, 6, ... , 2n} consists of all the even 
numbers. Let G be the union of all {d2j + 1) x (d2j + 1) minors of Xj for 
j G [n] such that d2j < min(d2j-i) ^2^+1) together with the union of all the 
2x2 minors of Yj^i for < j < n and i G [(i2j+i]. Then G is a Grobner 
basis for Ihaa- 

Proof. Note that lH,A,d reducible with separator S = {2n\. Thus, by 
applying induction and liberal application of Theorem 13. 7| we need only de- 
termine Grobner bases for Ih',A'4' ™ special case of A' = {{1, 2}, {2, 3}} 
and H' = 2 for all triples d' = (di, ^2, d^). Then we can lift the polynomials 
to get a Grobner basis for lA,H,d- In this special case, we are considering 
the ideal that is the kernel of the ring homomorphism 

(j)d : K[qi^,i^ I ii G [rfi], 23 G [ds]] ^ ^[ai^i2, hih I ij € [dj\] 

i2=l 

Thus, LH',A'.d is the vanishing ideal of Sec'^'^^^(F^'^^^ x F'^-'^^), the variety 
of secant {d2 — l)-planes to the Segre embedding of P'^i^i x F'^'-^^^. If d2 > 
min(di,d3) then lH,A,d = (0) is the zero ideal. If d2 < mm{di,ds) then 
Ih,A4 is generated by the (^2 + 1) x (^2 + 1) minors of the matrix {qi^,i^). 
Lifting these determinantal polynomials to yields the minors of matrix 
Xj. The minors of the matrices Yj^i are the elements of Quad^ for each of 
the toric fiber products that are used in building up lH,A d- ^ 

3.4. Group-based models on phylogenetic trees. In this section, we 
show how the toric fiber product arises in the construction of phylogenetic 
invariants for group-based models on phylogenetic trees. The fact that these 
phylogenetic models are toric fiber products plays a significant role in ^Jj 
where phylogenetic invariants for the group based models were originally 
constructed. The toric fiber product also plays a prominent role in [HI where, 
in the case of trivalent trees with underlying group Z2, it is proven that these 
phylogenetic ideals are Gorenstein, their Hilbert polynomials are computed, 
and their deformations are studied. For the sake of simplicity of exposi- 
tion, we will describe the underlying models in the Fourier coordinates (as 
opposed to the probability coordinates) and we only address the case of 
models whose labeling function is the identity map. We refer the reader to 
|17j and JH] for descriptions of these models in the probability coordinates 
and the application of the discrete Fourier transform. Also |17j contains the 
full description of these models with arbitrary friendly labeling functions. 

Let T be a tree with n + 1 leaves. Label the leaves 1, 2, . . . , n -|- 1, let the 
root of T be at the leaf n + 1, and direct the edges of T away from the root. 
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Given an edge e, a leaf / is called an descendant of e if there is a directed 
path from e to I. Denote by de(e) the set of all descendants of the edge e. 
Wc assume that the tree T has the property that for every edge e the set of 
descendants dc(e) is an interval of integers, i.e. de(e) = + 1, . . . j — 1, j} 
for some i < j & [n]. This amounts to saying that T has a drawing in the 
plane so that the leaves of T lie on a circle in numerical order. 

Let G be a group. Wc will use additive notation for G although G might 
not be abelian. If 31, ... is a sequence of elements in G, we denote by gg 
the sum of the group elements gi such that i G de(e); that is, 

9e= 9i- 
jede(e) 

Since G may not be abelian we use the convention that the sum is always 
taken in increasing order of the indices i G de(e). Let 

K[q] = K[qg^,„g„ \ gi e G] and K[a] = K[a^^^ \e G E{T), h E G] 
and consider the ring homomorphism 

(t)G,T ■ K[q] ^ K[a] 

Q9i...9n^ n 4f 

eeE{T) 

Definition 3.9. The ideal Ig,t = ^GT{(f)G,T) is the ideal of the group-based 
phylogenetic model with group G and tree T. 

Note that since T is an acyclic directed graph, there is an induced partial 
order on the edges of T. Namely e < e' if there is a directed path from e' 
to e. Let r be a tree that contains an interior edge (an edge not incident to 
any leaf). Then e induces a decomposition of T as the composition *T~ 
where T~ is the subtree of T consisting of all edges e' G T with e' < e and 
consists of all edges e' G T such that e' ^ e. Thus T~ and overlap 
in the single edge e. We root T~ by the tail of e, and keep the root of r+ 
at the original root n + 1. 

Without loss of generality, we may assume that the nonroot leaves of T~ 
consist of {1,2,..., k}, and the nonroot leaves of r+ are {e, k+1, ... , n}. Let 
]K[g]+ and IEC[g]_ denote the ambient polynomial rings of and Iqj^-, 

respectively. 

Theorem 3.10. Let T he a tree with an interior edge e, and resulting de- 
composition T = Tg+ * T~ . For each variable in K[q], K[q]^ and K[q]-, 
let deg{qg) = eg^, the standard unit vector with label g^- Then 

Ig,t = Ig,t+ ^^IcTr 
with A linear independent. 

Proof. Clearly A is linearly independent since it consists of standard unit 
vectors. To prove that Ig,t is a toric fiber product, we use our standard 
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technique of modifying the parameterization. As usual, it does not hurt to 
square a variable a^g] everywhere it appears. Thus we have: 



9e' 

e'eE{T) 



n n «^ 



,(e') rr a^/) 

e'&E{T+) e'£E{Tr) 



Thus, Ig,t is a toric fiber product. □ 

This allows us to deduce the main result from |17j . 

Corollary 3.11. If T is a trivalent tree the generators of Ig,t can be ex- 
plicitly determined from the generators of Iq^Ki 3 where i^i,3 is the three leaf 
claw tree. 

Proof. If T is a trivalent tree, it can be successively be decomposed by the 
* operation until each component tree is a -ftTi^s. □ 
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